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' Progress in the study of hadronic final states in deep inelastic scattering as well as pp, photoproduction 

^ ' and e'^e~ annihilation, as presented at the DIS97 workshop, is reviewed. 

CD ■ 
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X Introduction 

The large centre of mass energies and increasing statistical precision available at HERA, the Tevatron and 
LEP combined with recent theoretical developments open a new testing ground for QCD. The presentations on 
hadronic final states presented at the DIS97 workshop place increasingly significant experimental and theoretical 
constraints on the strong interaction. This summary focuses on the common themes and highlights of the 
hadronic final states working group. 



Fragmentation 



Perturbative QCD (pQCD) has proved to be a very successful theory in its application to hard processes. This 
enables the theory to be employed as a tool to tackle more complicated problems. One of these concerns the 
soft limit of QCD and colour confinement. Presently, multihadron production phenomena cannot be derived in 
a systematic way solely from perturbation theory without additional model-dependent assumptions. 

To investigate the limits of applicability of pQCD it is important to determine to what extent semi-soft phe- 
nomena in hard processes still reflect the properties of the perturbative evolution phase. This line of research, 
initiated almost 15 years ago has reached a high level of sophistication (see ^ for a recent review) and, as 
witnessed at this workshop, continues to inspire analyses in all major high-energy experiments. This summary 
deals with some of the most interesting results. 

Particle Rates 

In the current picture of hadron production, factorisation plays a predominant role in the different evolution 
stages of the process: 'preparation' of the primary partonic configuration, additional parton production de- 
scribed e.g. by (angular-ordered) parton showers (pQCD), hadron formation described e.g. by string or cluster 
fragmentation (non-perturbative QCD), secondary hadronic cascade-decays (QFD and non-perturbative QCD). 
To test such an ansatz, a comparative study of jet properties — including particle rates and spectra — in different 
reactions is required. At present this has not yet become a topic of primary interest at hadron-hadron colliders 
and HERA in spite of its importance for QCD. 

How much remains to be done is illustrated in Q where the impressive results from the LEP experiments are 
updated and reviewed. As far as hadron production is concerned, 38 different inclusive production rates of 
mesons and baryons are now measured at the and, for many of these, inclusive spectra are available. 

In all, good agreement is observed for the rates with tuned versions of JETSET 7.4 @ and HERWIG 5.9 %. 
A noteworthy exception is the baryon sector which remains an embarrassment for HERWIG. Either a better 
retuning or a critical re-evaluation of the cluster decay model seems required. The Lund JETSET string 
approach fares better but contains a large number of parameters related to flavour and spin. Since this number 
is increasing with time, little real predictive power is left. 

In e^e^ annihilation at LEP, evidence for breaking of 'jet universality' and factorisation may have been found 
from excess j] production — above JETSET expectations — at large momentum in three-jet events |^] (glueball 
production or surplus iso-singlet hadrons?) while no anomaly is seen for tt*^ production. It could be that 
the long-awaited direct manifestation of gluon jet fragmentation has finally been found If so, even larger 
discrepancies could be expected for the r/ @- The /o(975) and ao(980) mesons could also play a special role 
in the dynamics of quark confinement [[lO[ . A comparative study of these and other hadrons in quark and in 
gluon jets is called for. 

Problems also appear with strangeness production (mainly K and A) where DELPHI notes a deficit of strange 
particles in extreme two-jet events : the production of strangeness depends on the event topology in a manner 
that is not quantitatively described by JETSET. 

Although HERA experiments have only started to investigate the field so thoroughly explored at LEP, and 
information on identified particles is still scarce, first evidence has been found that the level of strangeness 
production in DIS and photoproduction, translated into a. s/u relative production rate is close to 0.2, to be 
compared with 0.3 in e+e^ 

Lessons to be learned from the vast amount of data in e+e" annihilations at the are that deviations from 
'universal fragmentation' may well have been observed and that the topology of the confining QCD fields is likely 
to play a role in hadroproduction. The rich variety of such topological configurations possible in ep collisions 
poses a real challenge for the experimentalists. 
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Particle Spectra 



A striking prediction of the perturbative approach to QCD jet physics is the depletion of soft particle production 
and the resulting approximately Gaussian shape of the inclusive distribution in the variable ^ = log E jet /E for 
particles with energy in a jet of energy Ej^t — the famous "hump-back plateau" Due to the intrajet 
coherence of gluon radiation, not the softest partons but those with intermediate energies (E oc E^^^^'^'*) 
multiply most effectively in QCD cascades. 

The shapes of the measured particle energy spectra in e^e~ annihilation turn out to be surprisingly close, over 
the whole momentum range down to momenta of a few hundred MeV, to the perturbative predictions based on 
the Modified Leading Log Approximation (MLLA) These observations can be taken as evidence that the 
perturbative phase of the cascade development indeed leaves its imprint on the final state hadrons. This, in 
turn, suggests that the conversion of partons into hadrons occurs at a low virtuality scale (of the order of the 
hadron masses), independent of the scale of the primary process, and involves only low-momentum transfers. 
This Local Parton-Hadron Duality (LPHD) may be connected to pre-confinement properties of QCD which 
ensure that colour charges are compensated locally LPHD remains, however, a strong hypothesis that 

is supposed to be valid only in an inclusive and average sense. With LPHD, only two essential parameters 
are involved in the perturbative description: the effective QCD scale A and a (transverse momentum) cut-off 
parameter Qq, resulting in a highly constrained theoretical framework; non-perturbative effects are essentially 
reduced to normalisation constants. 

New data on charged particle spectra were presented at this workshop by HI ZEUS |^ and CDF The 
HERA experiments concentrate on the current fragmentation region in DIS and perform the analysis in the Breit 
frame, where the exchanged boson is completely spacelikc. The new data confirm with much increased statistical 
significance the features observed in e+e~: approximately Gaussian shape of the ^ spectra with peak-position 
and width increasing with Q as predicted in MLLA. Moreover, for sufficiently large Q, they demonstrate the 
expected equivalence of the current region with one hemisphere of an e'^e~ event. With increasing luminosity 
being accumulated at HERA, this work should be extended to include moment and cumulant analyses of the 
spectra for which detailed predictions exist 

Beautiful confirmation of the MLLA-I-LPHD approach has been presented by CDF . This experiment studies 
charged particle momentum distributions in subsamples of dijet events. For fixed dijet masses (hence fixed jet 
energy) in the range 83 < Mjj < 625 GeV, the ^ distribution of tracks, within cones of various opening angle Q 
(with respect to the jet axis), is studied (see Fig. 0(a)). As dijet mass x jet opening angle increases, the peak of 
the spectrum, ^o, shifts towards larger values of ^ in perfect agreement with MLLA predictions and e+e^ data, 
as shown in Fig. |l|(b). Similar analyses should be possible in DIS and photoproduction at HERA but have not 
yet been attempted. 

Although present data on charged particles appear to confirm strikingly the perturbative approach to soft 
hadronisation, the situation is less clear-cut when spectra of identified particles/resonances are examined. At 
LEP, the conclusion is unambiguous: the peak positions do not agree with the naively expected mass dependence. 
Here also, data on spectra of different hadron species and from different jets at the Tevatron and HERA would 
be most helpful. 



Limiting behaviour at low momenta 

The analytical perturbative approach allows one to predict the limit of the one-particle invariant density in 
QCD jets Edn/d^p = dn/dy d^px at very small momenta p or, equivalently, in the limit of vanishing rapidity 



and transverse momentum |17|. If the dual description of hadronic and partonic states is adequate down to 
very small momenta, a finite, energy- independent limit of the invariant hadronic density, /q, is expected. This 
is a direct consequence of the colour coherence in soft gluon branching. A possible rise of Iq with centre-of-mass 
energy would indicate that either coherence or the local duality (or both) break down. Since colour coherence 
is a general property of QCD as a gauge theory, it is the LPHD concept that is tested in measurements of the 
soft hadrons. 
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CDF PRELIMINARY 



CDF PRELIMINARY 



Mjj = 390 GeV/c' 
MLLA FIT 



• = 0.456 Q„=234±20 MeV Const = 0-54±0.08 
□ 9 = 0.361 Q„=221 ±20 MeV Const = 0.52±0.0a 

• = 0.280 Q„=2 1 5±70 MeV Const = 0.51 ±0.08 
o = 0.217 Q„=2I4±20 MeV Const = 0.49 ±0.08 
A = 0.168 QErr=21 5±20 MeV Const = 0.47±0.08 
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(b) Evolution of the peak position with Afjj©. 



Figure 1: Comparison of preliminary CDF inclusive momentum distributions with MLLA predictions and €^ e 
annihilation data. 



The e+e^ annihilation data on charged and identified particle inclusive spectra have been found to follow the 
MLLA prediction surprisingly well, also at low centre-of-mass energies. The invariant spectra at low momentum 
scale approximately (within 10%) between 1.6 and 140 GeV and agree with perturbative calculations which 
become very sensitive to the strong running of as at small scales []l8| . 

At this workshop, HI presented the first Breit frame measurements of the invariant energy spectra in DIS as a 
function of Q, (see Fig. 3 in [Q). For sufficiently high Q, the data show that the low-momentum limit in that 
region of phase space is essentially independent of Q and indeed similar to that in e+e^ annihilation. 



Scaling violations 

Whereas the preceding paragraphs dealt mainly with the semi-soft limit of the hadron spectra, promising 
results from HI and ZEUS were presented to this workshop on the scaling violations at larger values of the 
scaled momentum Xp = ^vlQ [ P^ , P^ - In QCD, scaling violations of the fragmentation function are expected, in 
full analogy with scaling violations of structure functions, due to increased gluon radiation. This leads to softer 
particle spectra with increasing energy. In principle, the scaling violations at large allow a measurement of 
(Xs and have been exploited for that purpose in e~ annihilation. Whereas different e~ experiments must 
be combined to cover a sufficient range in centre-of-mass energy, this can be accomplished for ep collisions in a 
single experiment. 

HI data provides evidence for violation of scaling in the current region of the Breit frame. The corresponding 
ZEUS data, shown in Fig. ^, have also been compared with the CYCLOPS NLO calculation ||l^ incorporating 
fragmentation functions taken from e+e" [|5j. The data at large show a weak dependence on the input 
proton parton densities but a clear sensitivity to as- With analysis of more, already existing, data it should 
become possible to use the complete NLO calculations to extract as in DIS. 
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ZEUS 1994 Preliminary 
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Figure 2: 1/adcr/dxp measurement as a function of Q^. The preliminary ZEUS data are compared to data from 
e~^e~ annihilation experiments and the CYCLOPS NLO calculation of the inclusive charged hadron spectra. 



Multiple Interactions 

At this workshop, CDF presented preliminary results on the measurement of events where two scattering 
processes occur in the same event For distinct processes A and B, the cross section for this "double parton 
(DP) scattering" is given by ctdp — (TA'^B/'^cff, where acS is a process independent parameter. 

Events were selected with a relatively low transverse energy (Et > 16 GeV) "photon" trigger in conjunction 
with three jets with Et > 5 GeV. The separation of DP events from the underlying QCD background is 
determined by studying variables sensitive to decorrelation effects. In particular, the azimuthal angle between 
the two best-balancing pairs ( "photon" -l-jet versus dijet) is approximately flat for the DP signal and enables a 
statistical separation of events. A new feature of this analysis is that events with displaced vertices, where the 
jets are reconstructed from separate origins, are used to evaluate aAcrs directly and hence reduce the theoretical 
uncertainties. This allows the first relatively precise determination of the effective cross section: 

cTeff = (14.5± l.Tiai) mb. 

No ^-dependence is observed, within the uncertainty of ~ 20%. Assuming a uniformly dense ball of partons 
and using the measured inelastic pp cross section, one expects CTeff =11 mb. The measurement represents a 
milestone in the study of multiple interactions and provides the first significant experimental constraint on such 
processes. 

Multiple interactions, where two or more partons interact in the same event, represent a considerable uncertainty 
in the analysis of photoproduction events at HERA. In particular, the extraction of the gluon content of the 
photon at relatively low Xj requires careful modelling of these interactions, since they can contribute up to 
50% of the cross section at the relatively low Et values {Et > 6 GeV) measured so far. The Tevatron result 
should aid in a realistic estimate of the uncertainties due to multiple interactions in the extraction of the gluon 
density of the photon at HERA. Similarly, such measurements improve background estimates to di-boson and 
boson-|-jet production at the Tevatron as well as the predictions of jet rates from multiple interactions at the 
LHC. 
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Event Shapes 



The measurement of event shape variables has been well-estabhshed in e+e~ annihilation experiments. An 
important point in the development of our understanding of QCD is to ensure that the measurement is well- 
defined theoretically at the required level of precision. In this case variables are chosen which are relatively 
insensitive to soft gluon emission and coUinear parton branching. A determination of as(/i) is therefore possible 
by comparison with NLO theory plus resummed series or NLLA calculations. At this workshop, impressive 
results from LEPII were presented which enabled a LEP average as{Mz) = 0.120 ± 0.005 to be extracted | |2l|] . 

Recent theoretical developments in the understanding of infrared renormalon contributions, which lead to 
divergences from the perturbative calculations, allow the first steps to be made towards a direct comparison of 
theory and data without invoking hadronisation models. These power corrections, with a characteristic 1/Q 
dependence, have been calculated for event shape variables and could also be calculated for differential jet 
rates. 

HI presented new measurements of the thrust, jet broadening and jet mass in DIS for momentum transfers, 
7 < Q < 100 GeV, in the current region of the Breit frame j23|. The mean values of the event shape data show 
similar trends to results from e"'"e~ annihilation experiments as a function of Q. In the DIS case, one advantage 
is that the event axis is determined by the direction of the virtual boson, whereas in e~^e~ annihilation the 
axis has to be determined from the final state hadrons using e.g. the thrust axis. The data, shown in Fig. ^, 
have been fitted to NLO theory plus the calculated power corrections of Dasgupta and Webber. The important 
conclusion is that the size of the power correction, characterised by the parameter dg, is consistent with a single 
value of do = 0.491 ± 0.003 (exp) I^q (theory) for three of the four event shape variables. In the case of 
jet broadening, the calculation of the power corrections is subject to large uncertainties: hence this particular 
variable does not satisfy the requirement of being theoretically well-defined and is not included in the global 
fit. The development of these power corrections is not only intrinsically important, but should also enable more 
precise extractions of as by constraining the hadronisation uncertainties more precisely. 



HI op data 

fit to HI i-p d^ta 

NLO qCD faloulatio 




Figure 3: Mean values of event shape variables as a function of Q, from HI. The values are for 1-thrust, 
calculated using (a) the thrust axis or (b) the photon axis, (c) the jet broadening and (d) the jet mass. The 
dotted line indicates the NLO calculation. The full line indicates the fit incorporating power corrections. 
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Jet Shapes 



The shape of the transverse energy distribution of particles within a jet produced in various interactions allows 
the primary parton source of the jet to be identified. In addition, the data provide strong constraints on the 
coherence properties of the showering partons and enable tests of the universality of the fragmentation process. 
In an analysis from the ZEUS Collaboration ||2^, the jet shapes measured in photoproduction and DIS were 
compared with those from e^e^ annihilation and pp experiments. Jets are measured using the cone algorithm 
with a cone radius of 1. The jet shape, ip{r), is defined as the average fraction of the jet's transverse energy that 
lies within an inner cone of radius r. The distributions shown in Fig. ^ are therefore integral plots with ip(r) = l 
at r=l, whose rate of fall-off measures how broad the jet is. The data shown are for minimum jet transverse 
energies around 40 GeV. It is observed that the DIS and e+e^ data contain ~ 70% of their transverse energy 
within a sub-cone radius of 0.2, consistent with well-coUimated quark jets. In contrast, the pp data jets are 
rather broad, with only ~ 50% of their transverse energy being contained within the same sub-cone radius, 
consistent with predominantly gluon jets in this Ex range. 
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ZEUS Preliminary 
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★ ZEUS DATA DIS NC 1 994 (e* p) 
■ OPAL DATA (e* e') 

A CDF DATA (p p) 

• DO DATA (p p) 



Figure 4: Comparison of jet shape measurements from ZEUS(DIS), OPAL(e+e ), CDF and D0 (pp). The 
jet energy ranges are 37 < Ei^^ < 45 GeV, 35 GeV< £;j'=*, 40 < £;^^* < 60 GeV and 45 < < 70 GeV, 
respectively. 



Photoproduction data (see Fig. 2 in [Q) were also studied as a function of pseudorapidity and transverse 
energy. The observed changes in jet shape were reproduced in models which incorporate both direct and resolved 
photon processes provided that the resolved processes include the multiple interactions discussed above. NLO 
calculations from Klasen and Kramer [ psf determine the jet shape only at the lowest non-trivial order. In 
order to describe the data, an i?scp parameter is introduced which determines when two partons are merged 
into a single jet. The jet shape distribution is well described by NLO calculations with an i?sGp parameter 
which increases with increasing rapidity in the proton direction, but which is in the range 1.3 < i?scp < 1-8. 
Differential distributions of the average transverse energy in intervals of cone radius will enable i?sop to be fitted 
and provide further constraints on the models. 



High-E'r Jet Results from the Tevatron 

The Tevatron Collider provides a unique opportunity to study the properties of hard interactions in pp collisions 
at short distances. The production of jets at large Et and its comparison with perturbative QCD calculations 
are of interest as they can serve as a test of the elementarity of the partons. 

The CDF and D0 collaborations have measured jet cross sections over ten orders of magnitude in (Pa/dErdr] 
up to Et — 500 GeV, half way to the kinematic limit. The challenge of measuring such a steeply falling 
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spectrum is the understanding of the energy cahbration of jets. The highest Et jets are not directly cahbrated, 
resulting in large uncertainties. In this kinematic region the NLO calculations are well understood at the 10- 
20% level. However, precise knowledge of the parton distribution functions in the proton is required before 
firm conclusions can be drawn from the comparison of data and theory. Collider data can constrain the parton 
distribution functions in the proton and especially the gluon distribution at moderate x. Kosower presented 
a formalism to make such an extraction possible using NLO calculations, while minimising the amount of 
numerical computation involved | |2^ . 

The preliminary (published) inclusive jet cross sections as measured by CDF using the 1994-1995 (1992-1993) 
data sample in the pseudorapidity region of 0.1 < |?/| < 0.7 are compared to the NLO QCD in Fig. ^(a) [27| , |28|| . 
The latter are based on calculations by EKS with CTEQ3M |^ parton densities, renormalisation and 
factorisation scales /i = Ei^*' /2, and the standard Snowmass jet cone algorithm. The data and the prediction 
are in excellent agreement for Et < 250 GeV; at higher Et, however, the data lie significantly above the 
predictions. 
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(a) CDF data vs. theory 



(b) D0 data vs. theory 



Figure 5: Ratio between experiment and theory for the inclusive jet cross section as measured by CDF and D0. 



D0 presented updated inclusive jet cross sections in the region of |?7| < 0.5 with significantly reduced systematic 
uncertainties (by about a factor of two), coming from a re-evaluation of the jet energy scale corrections |2^. As 
shown in Fig. g(b), these results are in excellent agreement with NLO QCD over the entire Et range. D0 com- 
pares the data to NLO QCD predictions using JETRAD js^, the CTEQ3M parton densities, the renormalisation 
and factorisation scales /i = E™'^^ /2, and a modified Snowmass jet cone algorithm with i?sop=1.3. 

It should be noted that the D0 and CDF data have been compared to NLO QCD with slightly different input 
parameters which can introduce an i?T-dependent variation of ~ 10% on the theoretical predictions Also 
the two measurements probe different 77 regions. 

In order to compare directly the results of the two experiments, D0 also performed the analysis in the CDF rj 
region. Figure ^ shows the CDF data points as compared to a fit of the D0 data in the 0.1 < |r/| < 0.7 region. 
The error band corresponds to the D0 systematic error which is mainly due to the jet energy scale uncertainty. 
The CDF data lie above the D0 fit but are within the experimental uncertainties. For a more quantitative 
comparison between the two experiments, the correlations in the systematic uncertainties of the two data sets 
must be taken into account. 
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□ CDF 94-95 DATA (Preliminary) 




Band shows D0 systematic 



50 100 i50 200 



iOO 110 400 410 100 



Figure 6: Residual plot of the CDF data with a fit on the 0.1 < \ri\ < 0.7 D0 data. The band shown represents 
the D0 systematic uncertainty. 



The dijet angular distribution is an ideal tool to determine whether any possible excess of events in high-ET 
inclusive jet production is due to new physics effects. The angular distribution of the outgoing partons is 
strictly governed by the helicities of the partons participating in the hard process and is relatively insensitive 
to the parton densities. Any unusual contact interaction (with effective scale A) will flatten the centre of mass 
scattering angle distribution (or create an excess of events at low x)- The CDF published results on dijet 
angular distributions give a lower limit of A > 1.8 TeV. Figure ^ shows the recent D0 x distributions which 
are in good agreement with NLO QCD ||2^. Using these data, D0 rules out at 95% CL a model where quarks 
couple with a universal contact interaction of scale A ~ 2.1 TeV. 




Figure 7: Dijet angular distributions for D0 data compared to JETRAD for LO and NLO predictions with two 
different renormalisation/factorisation scales. 
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Dijet Production in DIS at HERA 



A major area of study at HERA is dijet production in DIS, which has earher been used to determine as and 
constrain the gluon density. Previously, data have been compared to NLO semi-analytic calculations using 
JADE-type algorithms. Now, with the flexible NLO Monte Carlo programs MEPJET H and DISENT [|| 
available, comparisons with various jet schemes are possible. Results were presented at this workshop using the 
cone ||,||, JADE !§ and kr ^ algorithms. 

In order to measure the cross sections, detailed comparisons with models incorporating parton showers/dipole 
chains and a hadronisation phase have been made. In general, these data are well described by the ARI- 
ADNE mi program and are reasonably well described by LEPTO or HERWIG je) . The next stage in the 
development of ARIADNE by the Lund group is the Linked Dipole Chain (LDC) model, which was reported 
at this wor kshop H]. 

A problem highlighted at the workshop relates to various attempts which have been made to correct to parton 
level in an attempt to determine the gluon density or the strong coupling constant directly. However, the 
relationship between the NLO partons and ARIADNE/LEPTO/HERWIG partons is far from clear and this 
introduces an uncertainty for theorists who wish to compare with published data. A presentation of the data 
corrected to hadron level is therefore required. 

The general observation in various analyses, with a range of different kinematic cuts, is that the measured dijet 
cross sections/rates tend to be higher than those predicted by the NLO calculations incorporating a default 
coupling constant and parton densities which describe the total DIS cross sectiond^ 




Figure 8: ZEUS preliminary dijet cross section as a function of ^, compared to LO and NLO predictions. 

This is illustrated in Fig. || where the ZEUS preliminary dijet cross section, corrected to parton level, is compared 
to the NLO and LO predictions The cross section is measured as a function of ^ = x{l + AIjj/Q^), the 
momentum fraction of the parton entering the hard scattering process. The data are ~ 30% higher than 
the NLO calculation and this difference persists after taking into account variations in calorimeter energy 
scale, jet energy resolution, the Monte Carlo used to correct to parton level, the input parton densities or 
the factorisation/renormalisation scale. However, the shape of the cross section is well described by the NLO 
calculations and this can be used to extract the power dependence of the gluon at low-^, ^g{£,) oc This 
results in a value of A = 0.38 ± 0.04 ± 0.18 at = 4 GeV^. 

A further development has been taken in the HI analysis. Using the fcr algorithm in the Breit frame a 
global fit has been performed of the HI and NMC DIS cross section measurements as well as the HI preliminary 
dijet rates. In order to account for hadronisation effects, an additional power correction term is incorporated 

^) A similar excess is observed in the ZEUS "resolved" photoproduction dijet data when compared to the NLO calculations for 
the lowest Et > 6 GeV and relatively low (0.3 < x°^^ < 0.7) data (see Fig. 2 in Q). 
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into the fit of the dijet rates. Ahhough the functional form of these power corrections has not yet been 
calculated, it is clear from the fits to the data that such a term is required. An empirical function h{x) = 
a + (3ln{x/Xo) +j\n^{x/xo) +dln'^{x/Xo) is introduced, where a,/3,7 and S are additional parameters in the fit 
and Xo — 10^'^. The additional contribution to the dijet cross section is determined as Aa{x,Q'^) = h{x)/Q'^. 
The fitted form of h{x) is shown in Fig. ^ together with the results of the global fit incorporating this power 
correction term. A calculation of the power corrections would therefore enable a simultaneous determination of 
its magnitude and provide further constraints on as as well as the parton densities. 
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Figure 9: Analysis of HI dijet data. 



BFKL-motivated measurements 



Recently, much interest has been focused on the small Bjorken-a; region, where one would like to distinguish 
BFKL from the more traditional DGLAP evolution equation One of the dominant Feynman graphs 
responsible for parton evolution in DIS is shown in Fig. The Xi denote the momentum fractions (relative 
to the incoming proton) of the incident virtual partons and pxi is the transverse momentum of emitted parton 
i. Such "ladder-type" diagrams with strong ordering in transverse momenta, ~ PtHY , but 

only soft ordering for the longitudinal fraction xi > X2 > ■ ■ ■ > Xn ^ x are the source of the leading log 
contributions which are summed in the DGLAP evolution equation p3[ . In the BFKL approximation, transverse 
momenta are no longer ordered along the ladder while there is a strong ordering in the fractional momentum 

Xn <C Xn-l < . . . < ~ Xjet- 

BFKL evolution can be enhanced and DGLAP evolution suppressed by studying DIS events which contain an 
identified jet of longitudinal momentum fraction Xjet = Pz{j) I Eproton (in the proton direction) which is large 
compared to Bjorken x |Q. Furthermore, tagging a forward jet with prij) — Q allows little room for DGLAP 
evolution while the condition Xjet S> x leaves BFKL evolution active. Assuming BFKL dynamics leads to an 
enhancement of the forward jet production cross section proportional to {xjet/x)'^'^^^, where is the BFKL 
pomeron intercept, compared to the 0{a^) QCD calculation with DGLAP evolution pGf . 
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Figure 10: Gluon ladder diagram contributing to jet production in DIS. The position and kinematics of the 
parton which can give rise to the forward jet is indicated. 



In Fig. |Tl|, recent data from HI and ZEUS | |4^ are compared with BFKL predictions |^ and fixed order 
QCD predictions as calculated with the MEPJET ]32|| program at NLO. The conditions prij) — Q and Xjet ^ x 
are satisfied in the two experiments by slightly different selection cuts. HI selects events with a forward jet of 
PtU) > 3.5 GeV (in the angular region 7° < 6{j) < 20°) with 

0.5 < Pt(j)VQ' < 2 , Xjet ^ Ejet/Eproton > 0.035 ; (1) 

while ZEUS triggers on somewhat harder jets of prii) > 5 GeV and ri{j) < 2.4 with 

0.5 < Pt(j)VQ' < 4 , Xjet = Pz{j)/Eproton > 0.035 . (2) 



Clearly, both experiments observe substantially more forward jet events than expected from NLO QCD. A very 
rough estimate of the uncertainty of the NLO calculation is provided by the two dotted lines, which correspond to 
variations by a factor 10 of the renormalisation and factorisation scales and A recent BFKL calculation 
(dashed lines) agrees better with the data, but here the overall normalisation is uncertain and the agreement 
may be fortuitous. Also, we recall that both experiments observe more centrally produced dijet events than 
predicted by the NLO QCD calculations. Whatever mechanism is responsible for the enhancement in central 
jet production may also play a role in the enhanced forward jet cross section. Clearly these issues must be 
resolved before the evidence for BFKL dynamics can be elevated to the status of discovery. 

The multiple gluon emission in ladder-type diagrams is also studied in jet-jet decorrelations at the Tevatron. 
DO presented preliminary results as a function of the pscudorapidity separation of the two leading jets in an 
event The measurement is compared to HERWIG and PYTHIA ||] parton-shower Monte Carlo simulations, 
and to BFKL predictions. The soft gluon emissions are expected to decorrelate the transverse energy (Et) 
and azimuthal angle (0) of the produced jets as the rapidity interval between them increases. HERWIG and 
PYTHIA simulations reproduce the observed decorrelation reasonably well. However, the leading-log BFKL 
resummation predicts a larger decorrelation while a NLO QCD calculation underestimates the decorrelation 
effects. Therefore, no clear conclusion on the question of BFKL dynamics can be drawn from the present 
Tevatron data. 
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Figure 11: Forward jet cross section at HERA as a function of Bjorken x within (a) the HI and (b) the 
ZEUS acceptance cuts. The BFKL result of Bartels et al. is shown as the dashed hnc. The sohd and 
dotted hues give the NLO MEPJET result for the scale choice = Hp = ^(0.5 kxf with ^ = 0.1, 1 and 10, 
which provides a measure for the uncertainty of the NLO prediction. 



Inst anions at HERA 



Perturbative QCD successfully describes hard scattering processes. Beyond these, nonperturbative processes 
are predicted by QCD as well, e.g. processes mediated via instanton configurations in the path integral js^ . 
Of particular interest in ep collisions are instanton processes which simultaneously produce n / light qZqR pairs 
and hence violate chirality by AQs = 2nf units. Ordinarily these processes are exponentially suppressed, by a 
factor exp[— 47r/Q:s]. In conjunction with multiple gluon emission, however, this suppression is ameliorated by 
a factor exp[— 47rF(a;')/as(^)], where the so-called "holy-grail function" F{x') equals unity at x' — 1. F{x') is 
known to decrease with decreasing x' , to about 1/2 at x' ~ 0.2 but is not reliably calculable for small values of 
x' ||. 

The expected fraction, f^-^\ of instanton induced events, compared to generic DIS events at the same x and 
Q^, depends critically on the shape of the holy-grail function at small x' . Expectations range between lO^*' < 
f^^'' < 10~^ if F{x') approaches a constant below a;^j„ = 0.2 . . . 0.3 |Q. Because several qq pairs and gluons are 
produced isotropically, the striking signature of instanton-induced events would be very high particle multiplicity 
and high average transverse energy deposition over a large region of the available phase space. 

HI reported on a search for such events (ssj. The best limits are obtained from the non-observation of events 
with large charged particle multiplicities. For 80 < < 220 GeV, limits of Z^^' < 0.5% have been set [53|. 
These limits are still about one order of magnitude larger than expectations from instanton calculations 
but they begin to probe the interesting parameter range. 
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VF+ Jets Production 



Hadronic production of W and Z bosons provides a clean probe of perturbative QCD calculations. Three anal- 
yses were presented from the D0 and CDF collaborations utilising the large sample of pp collisions accumulated 
from the 1994-1995 Tevatron Colhder run [El. 



The first analysis from the CDF collaboration was a measurement of the W/Z+ > n Jets cross sections for 
n = 1 — 4. Figure ^2|(a) shows the inclusive associated jet multiplicity distribution for W and Z bosons. 
The uncorrected data are compared to the CDF detector simulation incorporating the VECBOS LO QCD 
calculation plus HERWIG parton shower and hadronisation. The CTEQ3M |Q parton densities were used. 
The band in the theoretical predictions represents the effect of varying the renormalisation and factorisation 
scales from = -I- of the boson to the (pt)^ of the partons. Using the hard scale, JVP + p^, the LO 
QCD predictions are about a factor 1.7 lower than the data, for all jet multiplicities. On the other hand the 
predictions for the softer scale, (pt)^, are in better agreement with the data. 



The D0 collaboration reported results on the ratio of the production cross sections for W 
Jets, 7?.^'^, as a function of the minimum jet transverse energy, shown in Fig 



12(b). 



1 Jet to W -I- 
The data between 20 and 



60 GeV are consistently higher than the DYRAD NLO predictions by about a factor of two. This is a rather 
curious result since it is in a domain where one generally expects QCD to work well. 
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Figure 12: CDF and D0 results for W/Z+ Jets production. 



D0 also investigated color coherence effects in 14^-1- Jets events. For this study events with a W boson and 
opposing jet were selected and the distribution of soft particles around the colorless W boson and the jet (colored 
parton) was measured. The color coherence signal is observed by comparing the multiplicity distributions of 
calorimeter towers with Et > 250 MeV around the W and around the jet. It is concluded that both angular 
ordering and string fragmentation are needed in PYTHIA |5| to describe the data. 
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Conclusions 



The development of NLO calculations for a wide range of hadronic final state variables and the increasingly 
precise data from HERA, the Tevatron and LEP has provided a detailed testing ground for the strong interaction. 
At the DIS97 workshop, beautiful agreement of the data from LEPII with QCD was presented. At HERA and 
the Tevatron various chinks in the armour of QCD were identified. The detailed comparison of these data 
with the latest developments in the theoretical framework will determine whether the established paradigms 
are sufficient to understand the many facets of hadronic final state production reported at this workshop. 
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